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We develop a chiral SU(3) symmetric relativistic mean field (RMF) model with a logarithmic 
potential of scalar condensates. Experimental and empirical data of symmetric nuclear matter 
saturation properties, bulk properties of normal nuclei, and separation energies of single- and double- 
A hypernuclei are well explained. The nuclear matter equation of state (EOS) is found to be softened 
by tr^" mixing which comes from determinant interaction. The neutron star matter EOS is further 
softened by A hyperons. 
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I. INTRODUCTION 

The equation of state (EOS) of the dense hadronic 
matter is one of the keys in nuclear physics as well as 
in physics of compact stars In dense matter, var- 

ious forms of matter are expected to appear such as the 
hyperon admixture [3-[l0| , meson condensation [ll|, [l^l , 
baryon rich quark gluon plasma [isj , and color supercon- 
ductor Among these exotic forms of matter, hyper- 
ons are expected to appear at relatively small densities, 

(2 — 3)po- In addition to the variety of particle species, 
partial restoration of chiral symmetry is also expected in 
nuclear medium, and it vifould modify the properties of 
dense matter significantly. Therefore, it is desired to re- 
spect both hypernuclear physics information [l5l - [l9| and 
chiral symmetry (20| - [2^ in constructing the dense matter 
EOS. 

It is widely believed that hyperons should emerge such 
as in the neutron star core [4 [6l-[l0|. and/or during the 
black hole formation [2^ , as the baryon density increases 
in the hadronic matter. In neutron stars, hyperon ad- 
mixtures soften the EOS and reduce the maximum mass 
of neutron stars. They also increase the proton fraction, 
which may promote faster cooling processes [l^l . In black 
hole formation processes, hyperons are abundantly pro- 
duced due to temperature or density effects, and shorten 
the duration time of neutrino emission [23|. Hyperon 
admixtures are governed by the hyperon potentials in 
nuclear matter, which may be determined from the hy- 
peron separation energies from hypernuclei and hyperon 
production spectra. From this point of view, we need to 
adopt a theoretical framework which can explain both 
nuclear matter and finite nuclear properties. 

A chiral symmetry is another important ingredient in 
dense matter. It is a fundamental symmetry of QCD with 
massless quarks, and its spontaneous symmetry break- 
ing generates masses of constituent quarks and hadrons 
[20l |. Hadron properties and EOS would be modified in 
nuclear matter due to the partial restoration of the chiral 
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symmetry [21[. Thus it is preferable for hadronic many- 
body theories to possess the chiral symmetry, including 
its spontaneous breaking and partial restoration at finite 
densities. 

We have recently developed a chiral SU(2) symmet- 
ric relativistic mean field (RMF) model, abbreviated as 
SCL2 RMF, with a logarithmic a potential in the form 
of — log a |25| , which is derived from the strong coupling 
limit of lattice QCD (SCL-LQCD) If we naiVely 

include the vector meson in the linear a model ((f>'^ the- 
ory), the chiral symmetry is found to be restored below 
the normal nuclear density (Lee- Wick vacuum, or chiral 
collapse) [2^, [2^ . To avoid this problem, there are sev- 
eral attempts some of which result in having instability 
at large a values [13] or too stiff EOS [13. In the SCL2 
RMF, we do not have any instability in the a poten- 
tial, and the obtained nuclear matter EOS is found to 
be reasonably soft. In addition, the bulk properties of 
finite nuclei (binding energies and charge rms radii) are 
well described. Then it is desired to extend the SCL2 
to the SU/(3) version in order to describe hypernuclear 
systems. We expect that this extension enables us to get 
detailed information of various hypernuclei. 

In this paper, we introduce a chiral SU(3) symmet- 
ric RMF model, abbreviated as SCL3 RMF in the later 
discussion, as an extension of SCL2 RMF. Some of the 
model parameters are constrained by chiral symmetry 
through the hadron masses and vacuum condensates, and 
some of them are determined by the nuclear matter and 
finite nuclear properties. We also determine remaining 
parameters, meson- A coupling constants, by adopting the 
SU/(3) symmetric relation for vector couplings and fit- 
ting existing A hypernuclear data for scalar couplings. 
We show that we can reproduce the separation energies 
of single A hypernuclei (5'a) [isj and the A A bond en- 
ergy (ASaa) in aaHs [13 by choosing the coupling con- 
stants appropriately. The EOS of symmetric matter is 
found to be softened by the scalar meson with hidden 
strangeness, C, = ss, which couples with a through the 
determinant interaction representing the effects of Ua(1) 
anomaly [H], HI] . We also discuss the neutron star mat- 
ter EOS and neutron star maximum mass in the present 
SCL3 RMF. 
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It is generally preferable to derive the dense 
matter EOS from bare baryon-baryon interactions. 
Non-relativistic calculations based on the variational 
method [13, [s^ and the g-matrix (36l439j have been car- 
ried out based on realistic bare baryon-baryon inter- 
actions. These non-relativistic microscopic calculations 
only with two body forces do not reproduce the nu- 
clear matter saturation point, and three-body forces are 
found to be essential to explain the saturation property 
of symmetric nuclear matter [s^ - Ibq} . In Ref. [la Isof . 
the modern microscopic AW, YN, and YY interactions 
were examined with three-body forces in the framework 
of Briickner Hartree-Fock theory. The calculated results 
of the maximum neutron star mass suggest the impor- 
tance of the hyperon admixture and the three-body force 
in neutron star core. In the relativistic Briickner Hartree- 
Fock (RBHF) theory [10], empirical nuclear matter sat- 
uration is explained quantitatively. It should be noted 
that a part of the three-body force effects are taken into 
account in RBHF via Z-type diagrams which come from 
relativistic treatments, but RBHF results in Ref. [lO] do 
not include the bare three-body forces such as those via 
baryon resonances. The EOS in RBHF is approximately 
reproduced in RMF with non-linear uj interaction, which 
is introduced to simulate the high density behavior of the 
vector potential [4l|. In this work, we follow the latter 
stand point: We start from the Lagrangian with sev- 
eral model parameters and determine these parameters 
by fitting existing data. As a result, the scalar and vector 
potentials in a symmetric nuclear matter are found to be 
consistent with the RBHF results, then we expect that 
the results with hyperons would be also meaningful. 

This paper is organized as follows. In Sec. |TT1 we intro- 
duce a chiral SU(3) potential derived from SCL-LQCD 
in RMF as an extension from the chiral SU(2) potential. 
In Sec. mil we investigate the properties of symmetric 
nuclear matter, normal nuclei, and A hypernuclei, and 
fix the parameters in RMF model so as to reproduce em- 
pirical and experimental data. Then, we can anticipate 
neutron star matter EOS with an RMF model which can 
explain experimental data. Finally, we summarize our 
results and give an outlook, in Sec. IIVI 



II. CHIRAL SU(3) RMF MODEL 

An energy density as a function of the chiral conden- 
sate, abbreviated as a chiral potential here, describes the 
spontaneous chiral symmetry breaking and its partial 
restoration through the chiral condensates, and these chi- 
ral condensates determine hadron masses. In Ref. [25| . 
Tsubakihara and Ohnishi proposed to apply the loga- 
rithmic chiral SU(2) potential derived from the strong 
coupling limit of lattice QCD (SCL-LQCD) and devel- 
oped the SCL2 RMF model. In Subsec. MM we briefly 
summarize how to derive the chiral potential in SCL- 
LQCD [H, HI]. Extension to a chiral SU(3) potential is 
described in Subsec. Ill Bl and a chiral SU(3) symmetric 



RMF is introduced in Subsec. Ill CI 



A. Chiral SU(2) Potential from SCL-LQCD (SCL2) 

The lattice QCD action consists of the pure Yang-Mills 
part and the fermionic part. The pure Yang-Mills part 
is proportional to l/g^, where g is the bare QCD cou- 
pling. In SCL-LQCD {g —> oo), we can ignore the pure 
Yang-Mills action terms, and only those terms including 
fermions Sp are kept in the action [26j . The fermionic ac- 
tion with staggered fermions in the chiral limit is written 
in the lattice unit as, 

Sf[x,X,U] =^^Vt^{x) [x{x)Uf,{x)x{x + fi) 



-X{x + il)Ul{x)x{x)] , 



(1) 



where ?7^(x) = (^^ly-^o+xi-i i-a:,,-i represents the stag- 
gered factor. After integrating out link variables C/^ in 
the leading order of 1/d expansion [2^, we obtain the 
following partition function 



Z= j V[x,x.U]e^]y[-SF[x.X,U]) 



' j T^ix, X] cxp 



I ^c,f!{x)VM{x,y)M{y)f!a 



x.y.a^p 



= j I'lx,X,cr]exp(-S'<x[x,x,cr]) , 
5'^=^ X! cr(a;)a/3VA,/(x,y)cr(y)/3a 
+ Y <y{y)ai3VM{y,x)M{x)i3a 

x,y,a..l3 



(2) 



(3) 



The mesonic composites are defined as Aiapix) = 
X'f}{x)x'i,{x) , and the auxiliary fields Gap are related 
to the expectation values of the mesonic composites, 
{(Tap{x)) = —{Maf}{y))- In these equations, the su- 
perscript a denotes color and the subscripts a and 
/3 show the flavors of the quark fields. The lat- 
tice mesonic inverse propagator VM{x,y) is given as 
Vm{x, y) ^ Y^t, i^y,x+fi + Sy.x-fi) /4iVc- From the first to 
the second line in Eq. ([2]), the one-link integral formula, 
/ dUUabUlj^ = SadSbc/Nc has been used. 

Here, we substitute the auxiliary fields with the static 
and uniform scalar E^/j and pseudoscalar IIq,^ matrices 
as the mean field ansatz, 



O-apix) = -I- ie{x)Ilap 



(4) 



where e{x) = (_i)2:o-(-a;i-i-a;2-i-i:3^ Since fermions are de- 
coupled in each space-time point, the grassmann integral 
can be easily evaluated and the effective free energy is 
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obtained up to a constant as, 



ha 



i(tr [cjVmct]) 



N,{\og&ci{VMa)) 



= ^tr [mHi] - ^logdct [M^M] 



d 

''m'r 



(5) 
(6) 



where (• • •) denotes the space-time average, d = 4 is the 
space-time dimension, and M represents the meson ma- 
trix, M = 'SajS+iliap, in which the s{x) factor is replaced 
with unity in (Ta^. 

While the coefficients ba and are fixed in the lat- 
tice unit in Eq. ([5]), they depend on the lattice spacing 
and the scaling factor connecting the meson field and the 
quark condensate, which should be chosen for a and tt 
to be in the canonical form. Furthermore, n/ species of 
staggered fermions corresponds to Nf = ^rif fiavors, and 
the coefficient modification may not be trivial when we 
take Nf = 2 for SU(2) or iV/ = 2 + 1 for SU(3). Thus we 
stipulate them as parameters to obtain physical meson 
masses. 

In SU(2), meson matrix is given as M = (o'+«T-7r)/\/2. 
Requiring that the chiral potential has a minimum at 
cr = /tt and fitting the pion mass to^, one parameter rria 
is left as a free parameter. Then, the chiral potential is 
given as, 

= -y log (det MMt) -f ^ tr(MMt) - c^tj 



log(a2 



7r2) + 



2 



-2a./scL(7^) + 

J-IT 



2 2 



2 2 



/scL(a;) = log(l -x)+x + Y 



(7) 



(8) 



where ip^ = /tt — cr, and the explicit chiral symmetry 
breaking term — Ccrcr is introduced. We have omitted 
pion self-energy terms and constants in the third line in 
Eq. (O. Parameters a^, ba-, are given as. 



fTrm^ 



(9) 



With the above logarithmic cr potential, full chiral sym- 
metry restoration is suppressed because of the repulsive 
contribution of at small cr. The present treatment of 
SCL-LQCD is referred to as the zero temperature treat- 
ment, where diverges at cr —>• 0. In the finite tempera- 
ture treatment of SCL-LQCD 27|, the divergent behav- 



ior of disappears, while has a finite negative deriva- 
tive at cr 0. This finite negative derivative is enough to 
suppress full chiral restoration at finite density, since the 
nucleon Fermi integral contribution behaves as pscr^ and 
we always have a minimum at a finite cr value. Therefore, 
we suppose that the present chiral potential would be 
a good starting point to investigate cold nuclear matter 
and nuclei. 



B. Chiral SU(3) potential from SCL-LQCD (SCL3) 



In order to apply the logarithmic chiral potential to 
hypernuclear systems, it is necessary to include mesons 
with hidden strangeness (ss) such as and (f> in addition 
to mesons made of u and d quarks (cr, co and p). Here, 
we replace the meson matrix M with that of SU(3), 



Mil 



Oq + i7r+ + iK 



M 



- i-K^ M22_ + iK 

iK- R^ + iK° 
1 



+ i A'O 
C + iVs 



Mil = ^[icT + ir]) 
M22 = [(a + 11]) 



ITT 



")] , 
°)] ■ 



(10) 

(11) 

(12) 



In a similar way to the previous subsection, the chiral 
SU(3) potential in SCL-LQCD may be given as. 



n' y 

^ _ - — log (det M'M't) + - tr (MM^) 



CffCr - c^C + Vkmt 



(13) 



where the explicit chiral symmetry breaking effects are 
included as Co-cr and C(^C, terms. Since the strange quark 
mass is not small compared with f-^ and we have 
taken account of its effects also in the shift of the meson 
matrix, M' = A/ -f diag(0, 0, (5s). This shifted meson 
matrix plays the role of the constituent quark mass in 
Eq. Q and appears in the logarithmic term of the chiral 
potential. 

The Kobayashi-Maskawa-'t Hooft interaction term [32l . 
[3^ is represented in a form of the determinant of meson 
matrix, 



Vi 



KMT 



-d' (detM-HdetM^) 



(14) 



This KMT interaction represents the Ua(1) anomaly ef- 
fects. Without Vkmt, the above chiral potential is in- 
variant under Ul(3) x Ufl(3) transformation in the chi- 
ral limit, 5s = Ca = C(^ = 0. In the real world, Ua(1) 
symmetry is broken by the anomaly. Kobayashi and 
Maskawa [l^l proposed the above determinant interac- 
tion term, and this term is derived as the instanton in- 
duced quark interaction vertex by 't Hooft [s^ . 

Now we shall decompose the chiral potential Vy in 
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Eq. (|13p into meson mass terms and interaction terms. 
Fx = - a' logia'O + ^ + C') - d'a\ - c^a - c^C 

a 

1 22,1 22, T// ^ 



2/scL 



/SCL 




(16) 



where ^Pa = f-n — o' and <f(=fc~C show the deviation 
of cr and C from their vacuum expectation values, respec- 
tively, and Vcr^ denotes the interaction energy density. 
In the logarithmic potential, shifted vacuum expectation 
value of ( reads f^- = + 6s- The other meson fields 
than a and ^ are shown by (jja, and ma and (51^ represent 
their masses and interaction terms. We have ignored the 
third order term (</?(t)^</'C coming from the determinant 
interaction. This term does not change the chiral poten- 
tial significantly around the vacuum, but it makes the 
system unstable at large values of cr and ^. This is be- 
cause we do not have polynomial terms such as cr^ and 
C*, which stabilizes the chiral potential in the (p'^ theory. 
Compared with the case of SU(2), where all of tr(MAfl'), 
detM and detM'l' are proportional to the same combi- 
nation, + Tz'^ , we have several different terms from 
log(detM'Af't), tr(AfMl') and detM-hdetM^ in SU(3). 

In Eq. (|15p . the crC mixing appears in the quadratic 
form of the meson fields, thus we have to diagonalize the 
mass matrix to obtain observed cr and C meson masses as 



{fa fc) 



( ml 








ml ) 





Ml 




A/2 



fa 



(17) 



Five out of six (a', 6', c^, c^, d! and 5s) parameters in 
this chiral potential are fixed by fitting observed meson 
masses of m^, mK and Af^, and vacuum expectation val- 
ues of cr and C (/tt and /^). Relevant meson masses are 
related to the parameters (a', d! , 5s) as, 



ml 


=b' 


2a' 


2d'fc , 


(18) 


ml 


=b' 


2a' 

J TT 


2d'fc , 


(19) 


m]i 




V2a' 


- V2d'U , 


(20) 


ml 


=b' 


a' 

+ 772 ' 
•'C 




(21) 


u 




- 2d'U ■ 




(22) 
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FIG. 1: (Color online) Energy density in vacuum as a func- 
tion of a in the SCL3 model (solid curve) is compared with 
those in the linear a (0*, open-triangles), SCL2 (dashed 
curve), and TMl (dotted curve) models. 



We regard m^ as a model parameter, and give a' and 
6' — 2d' as in the case of the SCL2 model, 



_ J TT 



b' - 2d'fc =- 



(23) 
(24) 



We assign the observed ( as /o(980), and the parameters 
b', d' and 5s are determined to reproduce m^^ — 138 MeV, 
mx = 496 MeV and = 980 MeV. 

Coefficients of the linear terms in cr and C are deter- 
mined to reproduce the vacuum expectation values, 

=/.|6'-^-2d7c}=./x"^' , (25) 

cc =&7c . (26) 

Once we fix the parameters in the chiral SU(3) poten- 
tial, masses of other scalar and pseudoscalar mesons are 
determined as shown in Appendix [Xj Calculated masses 
of these mesons are tabulated in Table ID They are in rea- 
sonable agreement with experimental values except for k. 

In Fig. [TJ we show the energy density as a function of 
a. We compare the SCL3 results with those in SCL2 [25j . 
TMl [m and the linear cr {(f)'^) models. We adopt the 
parameter m^ — 690 MeV, which reproduces the bulk 
properties of normal nuclei as explained later, and opti- 
mal ipQ value is chosen for each cr. When we only consider 
the quadratic term in ip„, the energy density behaves as 
mlipl/2. Thus the energy density in SCL3 can be twice 
larger than the results in SCL2 and TMl, in which m^ 
is around 500 MeV, while the calculated results shows 
similar values around a = /t^- This is because the opti- 
mal ^ value is chosen and reduces the energy density in 
SCL3. In Fig. [21 we show the energy surface as a func- 
tion of ip^ and ip(^ . The optimal value of t/?^ is modified 



5 



TABLE I: Parameters and masses of ag, k, ry and ry' mesons as functions of m„. Parameters are determined by fitting Tr,K 
and C masses (m^ = 138 MeV, niK = 496 MeV, = 980 MeV), and vacuum condensate of a and ( (/^ = 92.4 MeV, 
/<; = 94.5 MeV). 



ma 




a'//^ b'/f^ d'/U Ss/U 


rrtag 






M^, 


Ma 


(MeV) 


(MeV) 








(MeV) 


(MeV) (MeV) (MeV) 


(MeV) 


630 


816.7 


1 1 OR 


/ O.U ( 




1108.4 


1000.8 


536.0 


1069.7 


321.7 


640 


818.3 


1 1 AA 


7*^ n^i 




1108.3 


1000.7 


535.7 


1062.4 


344.6 


650 


819.8 


11.81 


73.03 


23.06 0.4175 


1108.1 


1000.5 


535.4 


1054.9 


366.4 


660 


821.4 


12.20 


73.01 


22.67 0.4001 


1108.0 


1000.3 


535.0 


1047.2 


387.2 


670 


823.0 


12.59 


72.98 


22.28 0.3840 


1107.8 


1000.1 


534.6 


1039.4 


407.3 


680 


824.7 


12.98 


72.95 


21.88 0.3690 


1107.5 


999.8 


534.1 


1031.3 


426.7 
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826.3 


13.38 


72.92 


21.47 0.3550 


1107.3 


999.6 


533.6 


1023.1 


445.5 


700 


828.0 


13.79 


72.89 


21.06 0.3421 


1107.0 


999.3 


533.1 


1014.6 


463.8 


710 


829.6 


14.20 


72.85 


20.64 0.3300 


1106.8 
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532.5 


1006.0 


481.6 


720 


831.3 


14.62 


72.81 


20.21 0.3187 


1106.4 


998.6 


531.8 


997.1 


499.0 


730 


832.9 


15.05 


72.77 


19.77 0.3081 
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Exp. 




980±20 


672±40 


547.85 


957.78 


400-1200 




(Pa(MeV) 

FIG. 2: (Color online) Energy surface of chiral effective po- 
tential Vo-c at Pb = in SCL3. Points show chiral condensates 
[tpa and ifi^) at finite densities, Pb/po = 0, 0.5, • • • 5. 



from zero to a finite value by the ctC coupling from the 
KMT interaction, and the energy density is reduced by 
this mixing. 



C. SCL3 RMF model 

We incorporate the chiral SU(3) potential V^q dis- 
cussed in the previous subsections into the SU(3) RMF 
model. We consider the following SU(3) RMF La- 
grangian, which describes baryons which couple with a 
and C(= ss) scalar mesons, and w, p'^ (denoted by i?), 



and (j) vector mesons, 

i 

1 ni^ 1 m?: 

+ ^(c^.a;''f -Kc(^-,^C) . (27) 

M* - gaUPa - gciVC ) (28) 
C/f ^9u.^ uj'' + 9p^ R^ + ^-Y^eA'' , (29) 

where V^^{V = w, i?, 0) shows the field tensor of the vec- 
tor meson V . The uj^ term in Eq. (P7)) is introduced to 
simulate the behavior of the vector potential at high den- 
sities in the RBHF theory [13,1411. The term V^c,{lp„, ipi;) 
is the scalar self-interaction, and we adopt the chiral 
SU(3) potential V^q in Eq. (fTS]). 

In this Lagrangian, we have several model parameters 
to be fixed; in V^^, c^j for the lu self- interaction, and 
the meson-baryon coupling constants gmB- We assume 
that (1) nucleon mass is fully generated by the chiral con- 
densate, Mjv = Qrr N f-TT , (2) the vector couplings obey the 
SU/(3) relation [laSSli ^^-^ (3) nucleon does not couple 
with hidden strangeness mesons (C and (f>) psi. l43|. Sup- 
pressed N(j) coupling is understood in the Okubo-Zweig- 
lizuka (OZI) rule [43, while the scalar meson- nucleon 
coupling may violate the OZI rule [13] • Thus the 
assumption of g^jy ~ may be regarded as an working 
hypothesis. 
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TABLE II: Parameter set determined from saturation point 
of symmetric nuclear matter, binding energies and charge 
rms radii of normal nuclei, sa of single A hypernuclei and 
ABaa in AAHe. Input constants adopted in this paper are 
also shown. We adopt the saturation point, {po, Eo/A) — 
(0.150 im--\ -16.3 MeV). 



(MeV) Ss (MeV) g^N ffpjv g^A Qca 

690 32.81 11.95 294.9 4.54 3.40 5.17 




0.1 0.2 0.3 0.4 0.5 0.6 

Pb (fm"^) 



Under these assumptions, we have four meson- 
baryon coupling constants, guN,9pN,gaA,gcA, as 
model parameters. We have totally six parameters 
(mo- ,Cuj,gi^N, gpN , ffcrA 7 5cA ) ■ ^'^^ parameters relevant 
to normal nuclear properties {ma,Cui, guiN,gpN), first we 
give rricr and fix g^^N and Cuj by fitting the saturation 
point. Next from the binding energies of Sn and Pb iso- 
topes, and gpN arc obtained. The separation en- 
ergies of A hypernuclei mainly reflects the core-A po- 
tential depth, J7a ^ —30 MeV, thus the combination 
5ctA'Pct(po) + gcAVc(Po) is obtained from this fitting pro- 
cedure. Finally, the ratio of ga-A and g^A is determined 
from the AA bond energy in the double A hypernucleus. 
All parameters arc tabulated in Table HIl 

There are two points to be noted in view of the chi- 
ral symmetry in the above Lagrangian. We omit pseu- 
doscalar meson effects, and we do not require chiral 
symmetry in baryon-scalar meson couplings. In dis- 
cussing the in-medium modification of the chiral con- 
densates, the linear representation is more convenient, 
where scalar and pscudoscalar mesons appear as chiral 
partners. In the mean field treatment with parity fixed 
single particle baryon states, the expectation values of 
pseudoscalar mesons disappear, and we omit the explicit 
role of pseudoscalar mesons. Based on the chiral pertur- 
bation (ChPT) theory [45| and in the extended relativis- 
tic chiral mean field model [i^ , two pion exchange effects 
were examined in the relativistic nuclear energy density 
functional and they would modify the scalar and vec- 
tor coupling constants in a density dependent way. We 
here assume that two pion exchange effects arc repre- 
sented in the coupling constants and the self-interaction 
terms of scalar and vector mesons. In the linear rep- 
resentation, it is possible to construct an SU(3) chiral 
symmetric Yukawa coupling term of scalar and pseu- 
doscalar mesons with baryons [47| , but we have to assume 
baryons transform as nonet and only D-type coupling 
appears. For octet baryons having both D- and F-type 
Yukawa couplings, it is necessary to introduce two-types 
of baryons [48|, or to invoke the non- linear representa- 
tion [49l. [soj. which arc out of the scope in this paper. 



FIG. 3: (Color online) Energy per nucleon in symmetric 
nuclear matter as a function of the baryon density. Solid, 
dashed, dotted, dot-dashed, dot-dot-dashed curves show the 
results in SCL3(m<, = 690), SCL2 [H, TMl [H, NLl [H 
and NL3 [5^, respectively. Calculated results in variational 
calculation (FP) [U and RBHF [i^l are also presented by 
open squares and and open circles, respectively. 

III. FINITE NUCLEI AND NUCLEAR MATTER 
A. Nuclear matter and normal nuclei 

First we discuss the EOS of symmetric nuclear mat- 
ter. There are three relevant parameters, nia-, guiN 
and Cui- For a given rua- value, latter two are de- 
termined by fitting the saturation point (po, Eq/A) = 
(0.15 fm~^-16.3 MeV/A). In Fig. H wc show cal- 
culated energy per nucleon (E/A) in TMl, SCL2, and 
the present model (SCL3) with m„ = 690 MeV. When 
we adopt rria ~ 700MeV, the EOS in SCL3 is consid- 
erably softer than those in TMl and SCL2. We also 
find the SCL3 EOS is in good agreement with the vari- 
ational calculation results by Friedman and Pandhari- 
pande (FP) [s^l, especially at around pq. SCL3 EOS has 
rather soft incompressibility K ~ 211.0 MeV and this 
result is comparable with the empirical incompressibility 
K = 210 ±30 MeV [5l|. 

Nuclear matter EOS at several po has been probed in 
heavy- ion collisions js^-fsj] . In Fig. SI we show the region 
of pressures consistent with the experimental flow data 
analyzed by using the Boltzmann equation model [s^ . 
Danielewicz, Lacey and Lynch suggested the range of 
the incompressibility 167 MeV < K < 380 MeV in the 
density range 2po < Og < 5po (s^ . Other theoretical 
model calculations [SJ, HIl also explain flow data at AGS 
and SPS energies with K ~ 300 MeV. The calculated 
pressure in SCL3 is consistent with the pressure range 
suggested in Ref. (s^ . 

The EOS softening is caused by the C meson, which 
couples with a through the determinant interaction. In 
Fig. [21 we show the density dependence of the equi- 
librium point in {(pa-,f(^) plane. Equilibrium values at 
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FIG. 4: (Color online) Same as Fig. [3] but for the pressure. 
Shaded area shows the region of pressures consistent with 
the experimental flow data analyzed by using the Boltzmann 
equation model [s^ . 
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FIG. 5: (Color online) Same as Fig. [3] but for the scalar and 
vector potentials. 

Pb/ Po = 0^5 are shown with points. Wc find that the 
system evolves along with the valley where C, values are 
finite. This C, variation reduces the quadratic part of the 
mesonic energy density, while both ip^ and contribute 
repulsively in the higher order term, V^q. As a result, 
the energy gain from the scalar mesons is suppressed a 
little, and it leads to a smaller vector coupling to repro- 
duce the saturation point. Cancellation of smaller scalar 
and vector potentials leads to a softer EOS in SCL3. 

The density dependence of the scalar and vector poten- 
tials in SCL3 are found to be qualitatively consistent with 
the RBHF results HO] at low densities, < 0.3 fm"^. 
In Fig. [5l we show the scalar and vector potentials, 
U% = Mil - Mn = -gaNVcr and = gu^NOo in sym- 



metric nuclear matter, as functions of density in SCL3 
in comparison with those in SCL2, TMl, and RBHF. 
Scalar and vector potentials grow almost linearly with Pb 
at very low densities, and they are suppressed at higher 
densities in RBHF via the exchange and correlation, and 
the relativistic normalization In RMF models, the 
suppression is caused by the non-linear terms of cr and w. 
The RBHF results lie between SCL2 and SCL3, and the 
density dependence at < 0.3 fni^'^ is well reproduced 
with SCL3. Once these potentials arc given, similar EOSs 
are obtained at low densities; the difference appears from 
the non-linear terms, whose residual effects are small at 
low densities. At high densities {ps > 0.3 fm~'^), SCL3 
gives softer EOS than that in RBHF. This difference may 
come from the density dependence of the vector poten- 
tial. As shown in Fig. [SJ the vector potential in SCL3 
is suppressed more strongly than in RBHF. Since the 
term is introduced to mimic the density dependence in 
RBHF, it may be necessary to include other types of non- 
linear interaction terms for vector field provided that the 
density dependence in RBHF gives the convergent result 
in the hole-line expansion. 

For normal finite nuclei, binding energies per nucleon 
and charge rms radii are controlled by two parameters, 
TO(j and gpN- We determine them by fitting experimental 
data of binding energies and charge rms radii of some 
stable nuclei (^^C, ^^O, 4°Ca, ^^Ca, ^^Ni and ^OZr) and 
Sn and Pb isotopes. In Fig. \6\ wc show the calculated 
binding energies per nucleon (B/A) of C, O, Si, Ca, Ni, 
Zi, Sn and Pb isotopes, and B/A and charge rms radii for 
some stable nuclei are shown in Tables Hill and HVl where 
we also tabulate NLl NL3 [Hi, and non-chiral RMF 
(TM [23,|4l|) results. 

From these results, we find that SCL3 RMF model well 
describes the bulk properties of normal nuclei with the 
parameters shown in Table [III especially those of Sn and 
Pb isotopes which we can treat as spherical. Since B/A 
in these nuclei apparently reflects the character of EOS 
around po, our choice of the saturation point and adopted 
values of parameters (mcr = 690MeV, gpN = 4.54) is 
appropriate in explaining these data. 



B. A hyper nuclei 

From the very early stage [l,!!^,!!!], RMF models have 
been applied to A hypernuclei. Later, in Ref. [s^ [60j . 
tensor type couplings between vector meson and A were 
introduced so as to explain the small Is splitting. This 
problem was also examined from the view of baryon- 
mcson density dependent coupling [6]| . In addition to 
(T, Lj and p mesons, hidden strange meson fields (C and 0) 
were introduced in order to explain additional hyperon- 
hyperon interaction Q. In these works, hyperon-meson 
coupling constants are determined based on the flavor 
SU(3) or flavor-spin SU(6) symmetry, A hypernuclear 
single particle energies, hyperon potential depths in nu- 
clear matter, and suggestions from YN g-matrix. In 
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FIG. 6: (Color online) Binding energy per nucleon in normal 
nuclei from C to Pb isotopes. Solid, dashed and dotted curves 
show calculated results of SCL3, SCL2, and TM (TM2 for 
C to Ca isotopes and TMl for Ni to Pb isotopes) models, 
respectively. Experimental data are shown in open triangles. 



TABLE III: Experimental and theoretical binding energies 
and of stable nuclei. The results obtained from the SCL3 
model are compared with those obtained from TMl, [ZJ] 
TM2, [13] NLl, [13 NL3, [H and SCL2 [H models and 
with experimental data. 



BjA (MeV) 



Nucleus 


exp. 


SCL3 SCL2 TMl 


TM2 


NLl 


NL3 




7.68 


6.91 


7.09 




7.68 






16q 


7.98 


8.11 


8.06 




7.92 


7.95 


8.05 


2»Si 


8.45 


7.85 


8.02 




8.47 


8.25 




*«Ca 


8.55 


8.64 


8.57 


8.62 


8.48 


8.56 


8.55 


«Ca 


8.67 


8.58 


8.62 


8.65 


8.70 


8.60 


8.65 




8.73 


8.44 


8.54 


8.64 




8.70 


8.68 


«OZr 


8.71 


8.67 


8.69 


8.71 




8.71 


8.70 


"«Sn 


8.52 


8.51 


8.51 


8.53 




8.52 


8.51 


196pb 


7.87 


7.84 


7.87 


7.87 




7.89 




2°«Pb 


7.87 


7.87 


7.87 


7.87 




7.89 


7.88 



Refs. [57|,l58l, a; = gah/gaN ^ g^K/guN ^ adopted 
from discussions of tt, p and lo exchanges |57l or by fit- 
ting single particle levels of A hypernuclei [58|, while 
in Ref. Q, the ratio is set to be a; = 2/ 3 from light 
quark counting arguments. In Refs. 0, |60 ]. the vector 
meson-hyperon couplings are fixed from the SU(6) re- 
lation or the additive quark model results, e.g. g^jA = 
gujH ~ 2g^A ~ 2/3.gcjAf, while the scalar meson-hyperon 
couplings are determined from the A hypernuclear sin- 
gle particle energies or the hyperon potential depths in 
baryonic matter. In Ref. [s^, several sets of parameters 
are compared in the range guj^/ guiN = 0.18 — 0.63, and 
it is concluded that the A single particle energies are not 
enough to determine the meson- A coupling constants. 
Since the experimental information on Y N interaction 



TABLE IV: Same as Table Hn] but for charge rms radii. 



charge rms radius (fm) 



Nucleus 


exp. 


SCL3 SCL2 TMl 


TM2 


NLl 


NL3 


C 


2.46 


2.47 


2.43 


- 


2.39 


- 


- 


16 

o 


2.74 


2.63 


2.62 




2.67 


2.74 


2.73 


28Si 


3.09 


3.06 


3.04 




3.07 


3.03 




«Ca 


3.45 


3.43 


3.44 


3.44 


3.50 


3.48 


3.47 


^«Ca 


3.45 


3.46 


3.46 


3.45 


3.50 


3.44 


3.47 




3.77 


3.78 


3.77 


3.76 




3.73 


3.74 


9°Zr 


4.26 


4.26 


4.27 


4.27 




4.27 


4.29 


""Sn 


4.63 


4.61 


4.62 


4.61 




4.61 


4.61 


196pb 




5.48 


5.48 


5.47 




5.47 




20«Pb 


5.50 


5.54 


5.54 


5.53 




5.57 


5.58 



is limited, it is valuable to study the hypernuclear sys- 
tems with the RMF models including the chiral potential 
which can explain normal nuclear property. In previous 
studies [13,14^, chiral SU(3) symmetric RMF models are 
proposed. In the linear |47l | and non-linear [49j repre- 
sentation, various types of meson Lagrangian are com- 
pared. These models describe the normal nuclear prop- 
erties very well [i^ [s^l, while the hypernuclear prop- 
erties in these models arc not satisfactory from a phe- 
nomcnological point of view. Especially S and S hy- 
perons are considered to feel repulsive [62| and weakly 
attractive [63j potentials, but these features are not ex- 
plained yet. This may be suggesting the importance of 
other types of meson-hyperon couplings other than the 
Yukawa coupling or the SU/(3) breaking effects. 

In the present work, we study single- and doublc-A 
hypernuclei in the SCL3 RMF model, which has already 
shown to work well in normal nuclei and nuclear mat- 
ter as demonstrated in the previous subsection. We 
adopt the SU/(3) relation for vector coupling, and scalar 
meson-A Yukawa coupling constants arc chosen to fit the 
A separation energies of single A hypernuclei and the 
AA bond energy in ^aHc. In this treatment, the scalar 
meson-A coupling terms do not necessarily preserve the 
chiral SU(3) symmetry, but this may be an appropriate 
prescription at present because of the phenomenological 
problems in the chiral SU(3) RMF mentioned above. We 
omit the vector meson-hyperon tensor couplings, which 
does not affect the EOS in the mean field treatment. 

In single- and double-A hypernuclei, we have four ad- 
justable parameters, gah, gc,K, ffwA and g^A. Here, we 
assume that the vector couplings obey the lowest order 
SU/(3) symmetric relation jlSl |42||. 



C 



BV 



V2{gs tr {AQ tr {BB) + go tr {B {Af„, B}) 

+ gF tr(S[A/„,B])} 
V2{gs tr {My) tr (BB) + gi tr (BMyB) 

+ (?2tr(SBA/„)} , (30) 
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FIG. 7: (Color online) Separation energy of A. Parameters are 
determined by fitting Sa of a'C and }^C. Sa of p, d and / are 
wait-averaged considering the state numbers in each levels. 




-80 

0.2 0.4 0.6 0.8 

Pb (fm"^) 



FIG. 8: (Color online) Single particle potential of nucleon and 
A. Solid and dashed curves show the potential of nucleon and 
A, respectively. 



where go = {gi+g2)/'2, gp = (gi-52)/2, and the Lorentz 
indices and gamma matrices are assumed. In this form 
of vector coupling, g^^A and g^f^ are already fixed by the 
vector coupling constants with nucleons, g^jN and gpNi 
as 

guK = -qQuiN - -^gpN, 9<PA = [QuiN + S^pAf) . (31) 

Scalar coupling constants, g^A and g^^A, are then tuned to 
reproduce existing data of S'a and the AA bond energy, 
ABaa, observed in the Nagara event, p^ . 

In Fig. [71 we show the calculated results of S'a with the 
parameter set in Table HIl Here, we evaluate the zero- 
point kinetic energy, Ezpe, with a harmonic-oscillator 

r. 1/3 

wave function as Ezpe = j ■ 41^core MeV. The results 
of the 5a for p, d and / levels are the weight-averaged 
ones of the spin-orbit partners. 

The scalar potential of A is given in the form of linear 
combination of the coupling constants and chiral conden- 
sates, 

Ua{Pb) ^ -[gaAfaiPs) + 5CAVc(Pb)] • (32) 

In Fig. ISl we show the single particle potential for A, 
defined as the sum of the scalar and vector potentials, 
^^A = C^A + ^A' where the temporal component of the vec- 
tor potential [/° is defined in Eq. ([2^ . This sum roughly 
corresponds to the Schrodinger equivalent potential for 
A, U^^^ =Ul + {E/Ma)UI. As in the previous stud- 
ies |15l - [l7| . U\ amounts to be around —30 MeV at po- 
Calculated S'a values are very similar as far as the scalar 
potentials at pQ are the same. 

The Is splitting between p^^^ and in ^■^C is calcu- 
lated to be 900 keV in the present treatment. This result 
is larger than the observed small Is splitting [63|. A small 
Is splitting would be obtained by including the tensor 
type couplings between meson and A hyperon (59| , since 
this coupling directly corresponds to Is force when we 



translate the Dirac equation in RMF into a Schrodinger 
equivalent form. This coupling, however, reduces only 
the Is splitting and docs not change the average loca- 
tion of Is partners. Thus, we suppose that determined 
gaA and g^A are not affected dramatically. Two pion ex- 
change anti-Zs force is also suggested as the origin of the 
small Is splitting and examined by introducing density 
dependent coupling between meson and baryon (6]| and 
in the in- medium chiral SU(3) dynamics [65|. 

In Fig. ini we show the relation of g^A and ^ija de- 
termined by fitting Sa of and ]^C. This relation is 
roughly evaluated as g^A = 9 a A + 5c A /2 — 6. All the pa- 
rameter sets on the solid line can reproduce S'a of various 
single A hypernuclei similarly to the one shown in Fig. [T] 
We also examine Ai?AA in this parameter plane and find 
that the parameter sets in the gray shaded area explain 
the experimental A_Baa value of aa-Hs within the error. 
Combining these results, we obtain the set of coupling 
constants as {gaA, 5ca) ~ (3.40, 5.17) shown in Table [TTl 
which explains Sa and the central value of the experi- 
mental bond energy (A_Baa) simultaneously. 

The above set of coupling constants {go- A, Sca) deviates 
from a naive estimate, gaA/gaN = 2/3 and gaA/dcA = 
\pl. While the trA coupling is small, the scalar potential 
for A is around 2/3 of that for nucleons, and it is still 
dominated by a. The additional scalar potential comes 
from C, which does not couple with nucleons directly but 
appears from the aC, mixing generated by the KMT term 
in the Lagrangian. Since ipc^ evolves as (^^ ~ <^(t/2 in 
symmetric nuclear matter as seen in Fig. [2l the scalar 
potential of A behaves as ~ —gaA<f(T- The effective 
coupling 5cta is around 2/3 of go- at, then the smaller A 
scalar potential from a is compensated by the C, meson. 
The obtained gaA is close to 1/3 oig^Ni as suggested from 
the two pion exchange [57| . These observations may be 
suggesting the importance of pions and KMT term at 
finite densities. 
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FIG. 9: (Color online) and (7,;a suggested from experi- 
mental data. Solid line shows the relation between g^A and 
g(;A determined by fitting 5a of ^.^C and a^C, and gray shaded 
area shows the region which explains the experimental A_Baa 
value of AAHe within the error. The filled point represents the 
pair of gtrA and g^A, which explains Sa and the central value 
of the experimental bond energy (ABaa) simultaneously, and 
is adopted in the present work. 



C. Neutron star matter 

In previous subsections, we have fixed all parameters 
for nucleon and A by fitting the empirical saturation 
point of symmetric nuclear matter, and the binding en- 
ergies of normal nuclei and A hypernuclei. Now we apply 
the SCL3 RMF model to neutron star (NS) matter. 

In NS matter, we require neutrinoless /3-equilibrium 
and charge neutrality condition. 



(33) 
(34) 



B 



I 



The first equation relate the chemical potentials as, 
M)i = + Me = MA and fie = fJ-fj.- Under this equilib- 
rium condition, we calculate the energy density (e) and 
pressure (P), and we solve the Tolman-Oppenheimer- 
Volkoff(TOV) equation. 



dP _ [P{r) + e{r)] [M{r) + 4T:r^P{r)] 
dr ~ r{r - 2M{r)) 



(35) 



where M{r) denotes the mass inside the radius r. Here, 
we neglect nuclear crust for simplicity. 

In Fig. [TUl we show the NS matter EOS. We com- 
pare the results of SCL3 RMF model only with nucleons 
(SCL3) and with A hyperons (SCL3A). We also show 
the NS matter EOS in SCL2 [2|, TMl fU, NLl [El, 
NL3 [11] and lOTSY [I^, where A hyperons are not 
taken into account except for lOTSY. The lOTSY RMF 
model is based on TMl model and includes all octet hy- 
perons (A, E,S). We find that the NS matter EOS in 
SCL3 is softer than those in other RMF EOSs without 



hyperons. Especially, when we include A, A hyperons 
appear at around 2po and NS matter EOS including A 
hyperons becomes further softer. 

In Fig. [Til we show the results of neutron star mass 
as a function of the central density with SCL3, SCL3A, 
SCL2, TMl, lOTSY, NLl and NL3 models. While the 
EOS in SCL3 is much softer than other EOSs, calcu- 
lated maximum NS mass in SCL3 (without hyperon) is 
1.65M0, which exceeds the precisely observed NS mass, 
1.44 Mq [g^. When A hyperons are included, the maxi- 
mum mass is calculated to be 1.4OM0 with SCL3A and 
underestimates the observed mass. 

This underestimate is caused by the softer EOS, partic- 
ularly in the high pB region. It is suggested that extra re- 
pulsion coming from three-baryon interactions or string- 
junction model [67j which are repulsive for all baryon uni- 
versally are needed to surpass the known NS mass data 
in non-relativistic calculations [s^ 113, • In RMF mod- 
els, the EOS of nuclear matter is stiff enough, and extra 
repulsion is not generally required to support 1.44 Mq. 
In SCL3 and SCL3A, however, the incompressibility is as 
small as the empirical value and the pressure at high den- 
sity region is compatible with the estimate in heavy-ion 
collisions [52j . Thus, we encounter the same problem as 
in the non-relativistic calculations and have to consider 
additional repulsions which have a large effect at high 
densities. One of the candidates of this extra repulsion 
may come from the au coupling, such as in the term of 
(T^w^ [1^, [1^ . The (7UJ coupling results in reducing of the 
vector meson mass at high densities, and is found to give 
very stiff EOS when combined with the linear a model. 
Inclusion of such coupling may stiffen EOS in high pg 
region after re-fitting experimental data and it may solve 
the underestimation of the maximum mass of NS. 

Now we examine the effects of S hyperons in neutron 
star matter. In neutron star matter, it was believed that 
the S~ baryon would appear at the lowest density among 
the hyperons provided that the potential for S baryons in 
symmetric nuclear matter is similar to that for A [HI, @ . 
The strength of E-nucleus optical potential have been 
studied from the atomic shifts of S~ (toI . l7l| . which is 
sensitive to the attraction in -nucleus potential at nu- 
clear surface. In the inner region of nuclei, the analysis of 
the E~ quasi free production spectra from (tt", K'^) and 
{K~ , TT"*") reactions have yielded that repulsive E-nucleus 
potential is favored [g^- From these point of view, we 
employ the repulsive E potential [t^ which is suggested 
from SU/(3) relation, Eq. ([50)1 . Coupling constants, gj^s 
and g^s are given as 

1 \/2 

guT. = 2 iduiN + 9pn) 1 = — {9uN - 9pn) ■ (36) 

From the SU/(3) relation, should be equal to guT., 

(37) 



SU(3) ^ f, 



In order to explain the atomic shift data of E , how- 
ever, we need to adopt a smaller value of 5ps. In an 
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RMF model fit. Mares et al. showed that both of Si and 
Pb atomic shift data [zHl are well explained by adopting 
9pT.l9pN — 2/3 [llj. With the present Lagrangian, we 
have fitted the atomic shift data and have obtained 
the values = 1.97 {gpT./gpN = 0.434) and g^Y. = 3.16 
under the assumption of the naive quark counting cou- 
pling ratio for cr and C,^ ffo-s = V^g^s- The calcu- 
lated EOS, mass of neutron star, and particle fraction 

~ Pi/ Pb with S hyperons are shown in Figs. 1101 1111 
and in the middle panel of Fig. [TH respectively. In neu- 
tron star matter, E~ starts to emerge around ps ~ 0.4 
as a substitute of leptons because of its negative charge. 
Other E hyperons S"*") do not appear even at lOpo 
because of the repulsive potential and the negative charge 
chemical potential. The calculated EOS and maximum 
mass of neutron star arc not affected much. The repulsive 
potential suppresses the effects of S baryons compared to 
those of A, which plays a decisive role as a substitute of 
the dominant component, n. 

It is interesting to find that the starting density of 
hypcron (p^ ~ 0.4 fm~^) is much lower than those in 
previous studies 0413, which also adopt repulsive S po- 
tential in symmetric nuclear matter. In these works, E~ 
appears at much higher density. The main difference in 
the present work is the coupling strength with the p me- 
son. In Refs. d, [l3l, the SU(3) value {gps/gpN = 2) is 
adopted and the repulsive interaction from p is strong 
in high density neutron star matter, while the present 
coupling {gpT./gpN = 0.434) is much smaller and the re- 
pulsive potential from p is weaker. In order to demon- 
strate this point, we show the particle fraction results 
with gpT./gpN = 2 in the right panel of Fig. [TH We find 
that Tj~ appears only at ^ 1 fm~ , which is qualita- 
tively consistent with previous works [THlOj. 

One of the problems in the present SCL3 is that the 
c^j value (c^j ~ 294.9) is larger than those in TMl 
(c^ = 71.3075) and SCL2 (c^ = 200) models. The po- 
tential term of CijW'*/4 strongly suppresses w meson field 
especially at high p^ [4l|. The lower Cuj value an RMF 
model has, the higher neutron star maximum mass the 
EOS shows as seen in Figs. [TU] and [TT] When we re- 
duce this parameter by changing mo- value and re-fixing 
all the parameter in the way as we discussed, calculated 
maximum mass of neutron star should be 1.8 Mq on the 
condition of = 725MeV and c^^ = 75.66. With this 
choice, however, we cannot reproduce the binding ener- 
gies of Sn and Pb isotopes. Thus, in addition to chiral 
potential, the form and strength of vector meson poten- 
tial is also important and should be investigated further. 



IV. SUMMARY AND DISCUSSION 

In this paper, we have proposed a chiral SU(3) sym- 
metric RMF (SCL3 RMF) model, and examined its prop- 
erties in nuclear matter, normal nuclei, A hypernuclei 
and neutron star matter. We adopt a logarithmic chiral 
SU(3) potential, as the energy density as a function of 
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FIG. 10: (Color online) Energy per baryon in neutron star 
matter. Solid, short-dashed, bold-dashed, dotted, thin- 
dashed, dot-dashed, dot-dot-dashed curves show the results 
of SCL3, SCL3A, SCL2, TMl, lOTSY, NLl and NL3 results, 
respectively. 
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FIG. 11: (Color online) Same as Fig. [TO] but for the neutron 
star mass as a function of the central density. 



a at Pb =0, derived in the strong coupling limit of lat- 
tice QCD [13 . The Kobayashi-Maskawa-'t Hooft (KMT) 
determinant interaction term (32l . Issl ] is also introduced 
in order to take account of the Uyi(l) anomaly. Since 
the chiral symmetry relates the condensates and hadron 
masses, the number of parameters are reduced by intro- 
ducing this symmetry. After fitting tt, i^, fo(980) masses 
together with f-^ and /(;, we have only one free parame- 
ter, mcr, in the vacuum part. Under the assumptions that 
the nucleon mass are fully generated by the chiral conden- 
sate (Mjv = gcrN fir) and that the nucleon docs not couple 
with ss mesons, we determine four parameters relevant to 
normal nuclei {rricr, gi^N, gpN, Cc^), by fitting the empirical 
and experimental data of symmetric nuclear matter sat- 
uration point, binding energies and size of normal nuclei. 
For A hypernuclei, we assume that the vector couplings 
obey the lowest order SU/(3) symmetric relation [15ll43|. 
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FIG. 12: (Color online) Calculated particle fraction Yi of n, p, A, E , e , and fi . The other E hyperons (E" and E^) do not 
emerge in this baryon density range. 



and the remaining two parameters (gaAjScA) ^ire deter- 
mined by fitting the experimental data of the separation 
energies of single-A hypernuclei and the AA bond energy 
in the double-A hypcrnucleus, aa-Hc. 

We find that the SCL3 model well describes the sym- 
metric nuclear matter properties and the bulk properties 
of normal nuclei: The equation of state (EOS) is found to 
be softened by the (j( coupling generated by the KMT in- 
teraction, and the incompressibility of symmetric nuclear 
matter is found to be X ~ 210 MeV, which is consistent 
with the empirical value, K = 210 ± 30 MeV [EH. The 
EOS around po is i n ag reement with the results of vari- 
ational calculations |34| . and the pressure in the density 
region of 2po < Pb ^ 5po is in agreement with the esti- 
mates from heavy- ion collision data |52| . The density de- 
pendence of the vector potential is close to that in the rel- 
ativistic Briickner-Hartree-Fock (RBHF) calculation [i^l 
at low densities, Pb < 3po- The binding energies of nor- 
mal nuclei from C to Pb isotopes are reasonably well ex- 
plained except for the jj— closed shell nuclei. Single- and 
double-A hypernuclei are also well described. Separation 
energies of A in single-A hypernuclei, 5'a, are mainly de- 
termined by the potential depth, and the AA bond energy 
depends on the (A coupling ^.^a more strongly than on 
the a A coupling g^A- 

The calculated maximum neutron star mass in the 
present SCL3 RMF model underestimates the observed 
neutron star mass 1.44 A/0. This underestimate would 
originate from the soft EOS of nuclear matter at high 
densities. The vector potential is suppressed more 
strongly from the linear behavior, uj ^ gi^NpB/rn^, than 
the RBHF results, and the EOS of symmetric matter is 
softer than the results in the variational calculation [U 
and RBHF [13] at high densities. The suppression of the 
vector potential is caused by the w self-interaction term, 
—c^Lu^/A, whose coefficient is large in SCL3 compared 
with previous RMF models with this term. Since this 
term is introduced phcnomenologically to simulate the 



suppression of the vector potential in RBHF [4l| , it would 
be necessary to introduce other types of coupling, such 
as the scalar meson- vector meson coupling, a^uj"^ [28l[69j. 
The scalar- vector coupling acts to modify the in-medium 
vector meson mass |2ll . |72| |. It may be interesting to 
invoke the results of strong coupling lattice QCD with fi- 
nite coupling effects, where the plaquette contribution is 
found to generate the vector potential [T^I . Another pos- 
sibility is to introduce the repulsive three-baryon force, 
which is widely adopted in non-relativistic theories [3^ - 

US- 

In this paper, we examine how A and S hyperons af- 
fect the neutron star matter EOS based on experimental 
data. The isospin dependence of S potential in nuclear 
matter is found to be important for the composition at 
high densities. S hyperons are not included because the 
data are not enough to constrain the potential [6^ . Since 
S hyperons may further soften EOS at high densities, it 
is necessary to find the mechanism of re-stiffening at high 
densities in order to construct reliable and chiral SU(3) 
symmetric EOS including all these hyperons. Explicit 
role of pions [i^ [tI] is another important subject to 
study in terms of relativistic nuclear many-body prob- 
lems. The Zs-like potential [gl, [T^I from pion exchange 
would improve the binding energies of j j-closed shell nu- 
clei. In addition, tensor suppression may be play a criti- 
cal role in EOS at higher densities. 
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Appendix A: The masses of scalar and pseudoscalar 
mesons 



where we use same parameters, such as a', b' and d', 
defined in Sec. [Ill We have the mixing term of rj and 77^ 
mesons, thus one have to diagonalize their mass matrix 
to obtain vacuum masses. We tabulate calculated masses 
as functions of nia- in Table H] 



We show the formulae of the scalar and pseudoscalar 
meson masses except for cr, tt and K, which we have 
already shown in Sec. [Til When we adopt expectation 
value of each mesons and mass values of tt, K and fo as 
constraints, mass of oq, t] and r/' can be represented as a 
function of parameter m„ . These masses can be read as 



T]' 



6' 



J TT 

V2a' 



f-irf. 



- + V2d'U , 



c 



A-b'-^ + 2d'f^ 

J TV 

2 ,/ «' 
■'C 



to; 



2 



+ to2J + 



. 4£2 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 

(A6) 
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